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Abstract. We use the localized principle component analysis to detect deviations from scale 
invariance of the primordial power spectrum of curvature perturbations. With the technique 
we make uncorrelated estimates of the primordial power spectrum with five wavenumber bins. 
In the framework of a minimal ACDM model, using the latest cosmic microwave background 
data from the WMAP and ACT experiments we find that more than 95% of the preferred 
models are incompatible with the assumption of scale-invariance, but still compatible with a 
power-law primordial spectrum. We also forecast the sensitivity and constraints achievable 
by the Planck experiment by performing Monte Carlo studies on simulated data. Planck 
could significantly improve the constraints on the primordial power spectrum, especially at 
small scales by roughly a factor of 4. 
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1 Introduction 


Measurements of anisotropies in the cosmic microwave background (CMB) have played an 
essential role in constraining on basic cosmological parameters, especially in probing the dy- 
namics of inflationary phase in the early Universe [1]. The Wilkinson Microwave Anisotropy 
Probe (WMAP) satellite has measured the CMB over the full sky down to 0.2° resolu- 
tion [2, 3]. Measurements at higher resolution made with the Atacama Cosmology Telescope 
(ACT) [4] and the South Pole Telescope (SPT) [5] can provide us with complementary infor- 
mation about the early Universe on scales smaller than those probed by the WMAP satellite. 
The ACT experiment now measures fluctuations on scales from 0.4° to an arcminute. A 
combination of the WMAP data and the ACT data would improve the constraints on the 
cosmological parameters. 

Inflationary models with featureless potentials generically predict a primordial power 
spectrum of curvature perturbations close to scale invariant. Such models are usually pa- 
rameterized by an amplitude of spectrum A,, a spectral index n, and its running index a, 
as 


lIn Pr(k) = In As + (ns — Vin") 4a; Pes ; (1.1) 
ko ko 

where kg is a pivot scale. The parameterization is a Taylor expansion in the logarithmic 
amplitude and logarithmic wavenumber space around the pivot point. The special case with 
Ns = 1 and a, = 0 results in the Harrison-Zel’dovich (scale invariant) spectrum. In the 
slow-roll inflationary models, the spectral index and running index are first and second or- 
der in the slow-roll parameters respectively, and thus they are expected to be small. For a 
power-law parameterization (a, = 0), 99.5% of the preferred models are incompatible with 
the scale-invariant spectrum by using the 7-year WMAP data if tensor modes are ignored [3]. 
Even adding a running index, a slightly tilted power-law primordial power spectrum without 
tensor modes is still an excellent fit to the data [3]. Although, by combining the WMAP data 
with the ACT data, the running index prefers a negative value at 1.80, indicating enhanced 
damping at small scales, there is no statistically significant deviation from a power-law spec- 
trum [4]. Moreover, before claiming that a power-law spectrum is excluded, one should 
investigate extensions of the minimal ACDM model which could produce a similar effect in 
the CMB spectrum (but not necessarily in the large-scale-structure power spectrum). For in- 
stance, the marginal indication for enhanced damping in the small scale CMB spectrum could 
also be explained by extra relativistic degree of freedom [6]. Other simple extensions of the 
minimal ACDM model include small neutrino masses, a spatial curvature, a free primordial 
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helium fraction, etc. Here, we do not consider such alternatives, stick to the minimal ACDM 
paradigm, and investigate only the issue of the primordial spectrum beyond the power-law 
assumption. Indeed, motivated by theoretical models or features of the observed data, other 
various parameterizations of the primordial power spectrum have been considered: for exam- 
ple, a broken power spectrum [7] due perhaps to an interruption of the inflaton potential [8], 
a cutoff at large scales [9, 10] motivated by suppression of the lower multipoles in the CMB 
anisotropies [11, 12], and more complicated shapes of the spectrum caused by features in the 
inflaton potential [13]. 

Measuring deviations from scale invariance of the primordial power spectrum is a critical 
test of cosmological inflation. Either exact scale invariance or a strong deviation from scale 
invariance could falsify the idea of inflation. However, a strong theory prior on the form 
of the primordial power spectrum could lead to misinterpretation and biases in parameter 
determination. Some more general approaches have been proposed to reconstruct the shape 
of the primordial power spectrum from existing data, based on linear interpolation [14], 
cubic spline interpolation in log-log space [15], a minimally-parametric reconstruction [16], 
wavelet expansions [17], principle component analysis [18], and a direct reconstruction via 
deconvolution methods [19-21]. The first three approaches are sensitive to the overall shape 
of the spectrum while the last three reconstruct methods are sensitive to the local features 
in the spectrum. Therefore they are complementary and needed to cross-check each other. 

In this work we focus on the uncorrelated band-power estimates of the primordial power 
spectrum to measure deviations from scale invariance, based on the local principle component 
analysis introduced to study dark energy in [22]. We apply the method to the 7-year WMAP 
data [3] and in combination with small-scale CMB data from the ACT experiment [4]. In 
our analysis we adopt two main astrophysical priors on the Hubble constant (Ho) measured 
from the magnitude-redshift relation of 240 low-z Type Ia supernovae at z < 0.1 [23] and 
on the distance ratios of the comoving sound horizon to the angular diameter distances from 
the Baryon Acoustic Oscillation (BAO) in the distribution of galaxies [24]. Moreover, we 
generate mock data for the Planck experiment and then make forecast using the Monte 
Carlo simulation approach. As expected, Planck could significantly improve the constraints 
on the primordial power spectrum. 

This paper is organized as follows. In Section 2, we first describe the method and the 
data used in this analysis. We then apply the method to the seven-year WMAP data and in 
combination with the ACT data and present our results. In Section 3, using a Monte Carlo 
approach we analyze the sensitivity of the Planck experiment with respect to the primordial 
power spectrum. Section 4 is devoted to conclusions. 


2 Uncorrelated constraints from current observations 


We consider a spatially flat ACDM Universe described by the following cosmological param- 
eters 


{OQph?,Ah?, Os, T, A1, Aa, --, As}, (2.1) 


where Qh? and Q,h? are the physical baryon and cold dark matter densities relative to the 
critical density, h is the dimensionless Hubble parameter such that Ho = 100h kms~!Mpc7!, 
©, is the ratio of the sound horizon to the angular diameter distance at decoupling, and 7 
is the reionization optical depth. Since we do not consider extensions of the minimal flat 
ACDM model in this analysis, we fixed the primordial helium fraction and effective neutrino 
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number to their standard values, and did not introduce neutrino masses or a tensor modes. 
The primordial power spectrum parameters, A; = In [10° PR (ki)| (i = 1,2,.2;5), are the 
logarithmic values of the primordial power spectrum of curvature perturbations Pr(k) at 
five knots kj, equally spaced in logarithmic wavenumber between 0.0002 Mpc7! and 0.2 
Mpc7!. To reconstruct a smooth spectrum with continuous first and second derivatives with 
respect to Ink, we use a cubic spline interpolation to determine logarithmic values of the 
primordial power spectrum between these nodes. Outside of the wavenumber range we fix 
the slope of the primordial power spectrum at the boundaries since the CMB data place only 
weak constraints on them. Using ln Pr(k) instead of Pr(k) for splines ensures the positive 
definiteness of the primordial power spectrum at the expense of making the primordial power 
spectrum non-linear in the parameters. Otherwise, we must discard such steps in the Markov 
chain if the interpolating spline between the knots goes negative due to steep slopes [16]. As 
discussed in [15], the method is insensitive to local features in the primordial power spectrum, 
but is sensitive to the overall shape. 

The primordial spectrum parameters, A;, are correlated due to the geometric projection 
from the primordial power spectrum to the angular power spectrum and gravitational lensing. 
These correlations are encapsulated in the covariance matrix of the primordial spectrum 
parameters, 


© = (Ai — (Ai)) (Ay — (Aj), (2.2) 


which can be obtained by taking the average of the Markov chain and marginalizing over other 
cosmological parameters. The diagonal elements of the covariance matrix are the variances 
of A; and the non-diagonal elements represent corrections between the A; bins that slowly 
decrease with increasing bin separation. To eliminate these correlations, we employ the 
localized principle component analysis to construct a new basis, where the new parameters 
A; are uncorrelated [25]. This variant of the principal component analysis has recently been 
applied to probe the dynamics of dark energy [22, 26, 27]. We diagonalize the Fisher matrix 
F = C7}, so that 


F = ODO, (2.3) 


where O is an orthogonal matrix and D is the diagonalized inverse covariance of the trans- 
formed bins. The localized principle component analysis corresponds to the weight matrix 
W = OTD" 2O, which is usually normalized so that its rows sum up to unity. The weights 
are fairly localized in wavenumber since D!/2 is absorbed into O. With this choice, the un- 
correlated parameters can be now obtained by changing the basis through the weight matrix 
rotation, A = WA. When discussing our results, we will generally refer to these uncorrelated 
estimates. 

Data: We use the 7-year WMAP data (WMAP7) and in combination with the 148 
GHz ACT data during its 2008 season. For the WMAP data, we use the low-l and high-I 
temperature and polarization power spectra. We also consider the Sunyaev-Zel’dovich (SZ) 
effect, in which CMB photons scatter off hot electrons in clusters of galaxies. Given a SZ 
template it is described by a SZ template amplitude Asz as in the WMAP papers [2, 3]. 
For the ACT data, we focus on the band powers in the multiple range 1000 < / < 3000. 
Following Ref. [4] for computational efficiency the CMB is set to zero above | = 4000 where 
the contribution is subdominant, less than 5% of the total power. To use the ACT likelihood 
described in [4], aside from Asz there are two more secondary parameters, Ap and Ae. The 
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former is the total Poisson power at l = 3000 from radio and infrared point sources. The 
latter is the template amplitude of the clustered power from infrared point sources. We 
impose positivity priors on the three secondary parameters, use the SZ template and the 
clustered source template provided by the ACT likelihood package, and marginalize over 
these secondary parameters to account for SZ and point source contamination. We adopt 
two main astrophysical priors: the present-day Hubble constant Hy = 74.2 + 3.6 km s71 
Mpc! measured from the magnitude-redshift relation of 240 low-z Type Ia supernovae at 
z < 0.1 [23], and the distances ratios, r,;/Dy(z = 0.2) = 0.1905 + 0.0061 and r,/Dy(z = 
0.35) = 0.1097 + 0.0036, measured from the two-degree field galaxy redshift survey and the 
sloan digital sky survey data [24]. Here r, is the comoving sound horizon size at the baryon 
drag epoch and Dy is the effective distance measure for angular diameter distance. 
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Figure 1. Uncorrelated constraints on the primordial power spectrum of curvature perturbations 
and their weight functions, derived from the WMAP7+Ho+BAO combination (top panels) and from 
the WMAP7+ACT+H,)+BAO combination (bottom panels). The blue and red error bars show lo 
and 2o uncertainties respectively. 


Results: Our analysis is carried out using a modified version of the publicly available 
CosmoMC package, which explores the parameter space by means of Monte Carlo Markov 
Chains [28]. Figure 1 shows the uncorrelated constraints on the primordial power spectrum 
of curvature perturbations (68% and 95% CL) and the corresponding weight functions that 
describe transformation from correlated parameters A; to the uncorrelated A;, derived from 
WMAP7+Ho+BAO (top panels) and from WMAP7+ACT+H)+BAO (bottom panels), re- 
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spectively. We can see that the power spectrum is best determined around k ~ 0.007 Mpc7!, 


and less accurately determined at much lower and much higher wavenumber because of the 
cosmic variance and dominant noise respectively. As shown in the top-left panel of figure 1, 
95% of the preferred models are incompatible with the assumption of scale-invariance, but 
still compatible with a power-law primordial spectrum. Adding the ACT data we find that 
there is more deviation from a simple scale-invariant spectrum due to reduced errors and 
suppressed spectrum at high-k, but it is weaker than the corresponding result from WMAP 
adopting a inflation-motivated power-law spectrum prior [3]. Note that the weights are fairly 
localized in k, as found in the context of dark energy measurements [22, 26, 27]. Moreover, 
the weight functions in the top-right panel are similar to those in the bottom-right panel. 


Parameter WMAP7+Hp+BAO WMAP7+ACT+AH_)+BAO Planck 

Aj 3.0776 + 0.3973 3.0696 + 0.3855 3.1488 + 0.3212 
Apo 3.1302 + 0.0720 3.1388 + 0.0681 3.1289 + 0.0379 
Ag 3.2451 + 0.0356 3.2532 + 0.0393 3.1327 + 0.0243 
A, 3.0515 + 0.0584 3.0872 + 0.0473 3.1366 + 0.0122 
As 3.0237 + 0.1076 2.9390 + 0.0784 3.1384 + 0.0167 


Table 1. Uncorrelated constraints on the primordial power spectrum with 68% confidence levels. 


3 Planck forecast constraints 


In this section, we apply Monte Carlo Markov Chain methods to assess the accuracy with 
which the primordial power spectrum can be constrained from Planck experiment. Following 
the approach described in Ref. [29], we generate synthetic data for the Planck experiment and 
then perform a systematic analysis on the simulated data. Assuming a fiducial ACDM model 
with a scale-invariant power spectrum, one can use a Boltzmann code such as CAMB [30] to 
calculate the angular power spectra Cr x Cc ; QFE , Ge and cP d for the temperature, cross 
temperature-polarization, polarization, deflection field and cross temperature-defection. We 
assume that beam uncertainties are small and that uncertainties due to foreground removal 
are smaller than statistical errors. For an experiment with some known beam width and de- 
tectors sensitivity, the noise power spectrum N, cae N fe and Ne can be estimated. Here we 
use the FuturCMB package! to calculate N aa based on the quadratic estimator method pro- 
posed in [31], which provides an algorithm for estimating the noise spectrum of the deflection 
field from the observed CMB primary anisotropy and noise power spectra. For Planck we 
combine only the 100, 143 and 217 GHz HFI channels, with beam width 6pwum = (9.6’, 7.0’, 
4.6’) in arcminutes, temperature noise per pixel op = (8.2, 6.0, 13.1) in uK and polarization 
noise per pixel og = (13.1, 11.2, 24.5) in uK (see Ref. [32] for the instrumental specifications 
of Planck). Given the fiducial spectra C, and noise spectra N;, one can generate mock data 
C;. We perform a Monte Carlo analysis through the likelihood function defined as 


ee RPDA (aeS) , (3.1) 


1The FuturCMB package is available at: http://lpsc.in2p3.fr/perotto/ 
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where 


D = ĈITOPEGH + CPTCPEGH + CTTOPP Ôg 


-OFP (CEEC# + 2CFECH) — OF? (CPACPF + 2CF ACP) , (3.2) 
_ a te a 255 A 
IG] = OFTCPECH — (CT®)” Cf4 — (CF4) CFF, (3.3) 
a x x à x 2 de A 2 x 
IC| = CFTCPECH — (CTP) êr - (CP) êP”. (3.4) 


Here, C; = Cı + N; is the theoretical spectrum plus noise and fsky is the sky fraction due 
to foregrounds removal. For Planck we choose fsky = 0.65, corresponding to a +20° galactic 
cut, and consider the data up to l = 2000. 
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Figure 2. Uncorrelated constraints on the primordial power spectrum of curvature perturbations 
and their weight functions from Planck simulated data. The blue and red error bars show lo and 20 
uncertainties respectively. 


Our results are presented in figure 2 and table 1 for Planck simulated data. As we can 
see in table 1, Planck will reduce the uncertainties in A;, especially in Ay by a factor of 3.9 and 
in As by a factor of 4.7. Since large uncertainties of the power spectrum at low-k mainly arise 
from the cosmic variance, measurement of A, is limited. The weight functions in the right 
panel of figure 2 are a little better localized in wavenumber than those in figure 1 because 
the weak lensing effect is extracted from CMB maps provided by Planck. Furthermore, we 
have checked that the weight functions depend weakly on the fiducial cosmological model. 


4 Conclusions 


Most inflationary models predict small deviations from a scale-invariant power spectrum. 
Therefore, the measurements of deviations from an exact scale-invariant spectrum would 
provide a firm probe of the dynamics of an inflationary phase happened in the early Universe. 
The local principal component technique is a powerful tool for measuring deviations from the 
scale-invariant spectrum, complementary to other approaches to reconstruct the primordial 
power spectrum or the direct testing of slow-roll inflation [33-35]. In this paper, we have 
used the localized principal component analysis to produce uncorrelated estimates of the 
primordial power spectrum of curvature perturbations. In the framework of a minimal ACDM 
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model, we found that more than 95% of the preferred models are incompatible with the scale- 
invariant spectrum, but still compatible with a power-law primordial spectrum by using the 
7-year WMAP data in combination with the ACT data. This conclusion is a little stronger 
than the corresponding result in Ref. [15], but weaker than when the inflation-motivated 
power-law prior is adopted. We have performed a systematic analysis of the future constraints 
on the primordial power spectrum achievable from the Planck experiment. We found that 
Planck would be able to shrink the error bars on the spectrum bins especially at small scales 
by roughly a factor of 4, which is promising to definitively detect these deviations. 
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